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In the first year of the grant, potential energy surfaces (PES’s) have been calculated 
for a number of reactions. The reactions currently under study are relevant to the de- 
sign of the SCRAM jet engine for the NASP. Two hydrogen/oxygen reactions have been 
studied. The first reaction, 


H + 0 2 -> H0 2 -*0H + 0 (1) 

is the most important chain propagation reaction in combustion. Although there have 
been previous studies of reaction (1), they focused on the H + O 2 region of the surface 
and a detailed examination of the OH + 0 region was needed. This latter region of the 
surface is complex, because chemical bonding and long range electrostatic forces com- 
pete, leading to an OH - 0 minimum at large 00 distances (% 6.0 ao) followed by a 
saddle point before the H00 well. Approximately 150 points have been obtained which 
characterize the minimum energy path(MEP) for the entire PES. Additional points will 
be computed to characterize the hydrogen atom transfer process. In collaboration with 
Ron Duchovic (Eloret Institute), the complete set of points will be used to generate a 
global PES suitable for dynamical studies. A more detailed discussion of this work is 
given in ref. 2. 

An additional hydrogen/oxygen reaction which has been studied is the reaction 

0( 1 D) + H 2 -» H 2 0* -» OH + H. (2) 

The reason for studying this system is that 0( J D) formed by photolysis of O 2 could be 
used as a means of initiating combustion at lower temperatures in the SCRAM jet en- 
gine. As discussed in the semiannual report for the first six months of this grant, both 
the vibrational energy distribution in the OH product and the OD/OH ratio when 0( 1 D) 
is reacted with HD seem to depend critically on the long range behavior of the potential. 
The new, more accurate, potential energy surface may explain discrepancies between 


previous dynamical studies and experiment. A more detailed discussion of this work is 
given in ref. 3 . 

Another important group of reactions for design of the SCRAM jet are three body 
recombination reactions. As a start in this area, calculations were carried out for the 
reactions 


H + N 2 -+ HN 2 ( 3 ) 

and 

H + HN 2 -4 H 2 + N 2 ( 4 ) 

since the combination of eqns ( 3 )-( 4 ) constitutes a mechanism for H atom recombination 
with N2 acting as a third body. HN2 is also a very important intermediate in Thermal 
DeNO* processes. Specifically, addition of NH3 leads to substantial NO* reduction in 
hydrocarbon combustion. The key step here is thought to be 

N E 2 + NO — ► products 

where the expected products are HN2 4- OH and N2 + H2O. The observed branching ra- 
tio varies from study to study and while there is no direct evidence that HN2 is formed, 
OH is observed but H atoms are not. The presence of an HN2 intermediate is also im- 
plicated indirectly by modeling studies which require chain propagation steps producing 
OH in order to explain the DeNO* process. 

A key factor, currently undetermined, is the lifetime of the HN2 species. The potential 
for dissociation of HN2 has a shallow well « 3 kcal/mole above the H -f N2 asymptote 
and separated from that asymptote by a barrier. The dissociation process is primarily 
controlled by quantum mechanical tunneling. Enough points on the HN2 dissociation 
potential have been obtained to characterize the MEP. Approximate calculations (one- 
dimensional tunneling through an Eckart barrier) have been carried out to determine 



the rate of dissociation via tunneling. Additional electronic structure calculations with 
larger basis sets were carried out in collaboration with Celeste Rohlfing( Sandia National 
Laboratories, Livermore) to assess the accuracy of the computed potential. The larger 
basis set calculations indicate that the PES is converged to « ± 1.0 kcal/mol. Based on 
this and the tunneling calculations the lifetime of HN 2 is estimated to be less than 10~ 8 
sec. This is much shorter than values of about 10~ 4 sec which would fit current models 
of thermal de-NO* processes. This result suggests that there may be other important 
processes not included in the current models. A more detailed discussion of this work 
will be found in ref. 4. 

Part of the potential for reaction (4) is shown in Fig. 1. Here the geometry of the 
HN 2 molecule is fixed at close to the computed equilibrium geometry and the distance 
between the second H and the center of mass of the N 2 as well as the polar angle are 
varied. From Fig. 1 it is seen that there are three channels which lead to stable minima 
( cis and trans HNNH and H 2 NN) and a fourth abstraction channel which leads to N 2 
+ H 2 . These results indicate that HN 2 is a very reactive species which is consistent with 
the short computed lifetime. At the current time it appears that the lifetime of HN 2 is 
too short for H atom abstraction from this species to contribute significantly to the H 
-f H + N 2 three body recombination process; i.e. the dominant process may involve a 
different mechanism such as stabilization of a metastable H 2 species. Study of the lat- 
ter process would require regions of the surface which were not considered in the present 
study; although, the HN 2 potential discussed above would be a good starting point. 

NO will be an important species in the SCRAM jet engine. Since H atoms will also be 
present the potential surfaces for HNO are important. Currently studies are underway 
to characterize the lowest three surfaces ( ] A', 3 A", and 1 A") of HNO. Fig. 2 shows the 
potential for H atom addition to the N end of NO. The lowest surface, which involves 
addition to the radical orbital of NO, has no barrier, while the other two surfaces, which 
involve addition to the 7r bond of NO, have barriers. The HNO region of these three sur- 
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Fig. 1. The potential for the interaction of an H atom with HN 2 . The HN 2 geometry is 
fixed at near the computed equilibrium geometry and the energy is shown as a function 
of i)the distance between H and the center of mass of N 2 and ii)the polar angle (0). 
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faces has been characterized. The computed geometries for the ground state and 
1 A" excited state as well as the computed excitation energy are in good agreement with 
experiment. Calculations are also underway to examine the NOH region of the same 
three surfaces. These results will be combined with calculations in collaboration with 




Celeste Rohlfing and Carl Melius at Sandia National Laboratory, Livermore on the 0 + 
NH and N + OH regions of the potential to produce a global potential for HNO. 

Fig. 2. The energy as a function of tjvh along the minimum energy path for H atom 
addition to the N end of NO for the lowest 1 A', 3 A", and 1 A" surfaces of HNO. 



In the next year a number of projects will be completed. As mentioned above, addi- 
tional points will be obtained for reaction (1). Studies on HNO will also be completed 
in collaboration with scientists at Sandia National Laboratory, Livermore. In collabora- 
tion with Ron Duchovic (Eloret), global representations of these potential surfaces will 
be constructed, and dynamical studies will be carried out using both statistical methods 
(RRKM theory and transition state theory) and classical trajectories. Further, quan- 
tum dynamical studies will be completed in collaboration with Dave Schwenke (Eloret). 
Finally, new studies using the improved 0( 1 D) + H 2 potential energy surface will be 
carried out by George Schatz (Northwestern University). 
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Abstract 

The minimum energy path (MEP) for the addition of a hydrogen atom to N 2 is 
characterized in CASSCF/CCI calculations using the [4s3p2dlf/3s2pld] basis set, 
with additional single point calculations at the stationary points of the potential 
energy surface using the [5s4p3d2f/4s3p2d] basis set. These calculations represent 
the most extensive set of ab initio calculations completed to date, yielding a zero- 
point corrected barrier for HN 2 dissociation of ~ 8.5 kcal mol -1 . The lifetime of 
the HN 2 species is estimated from the calculated geometries and energetics using 
both conventional Transition State Theory and a method which utilizes an Eckart 
barrier to compute one-dimensional quantum mechanical tunneling effects. This 
study concludes that the lifetime of the HN 2 species is very short, greatly limiting 
its role in both termolecular recombination reactions and combustion processes. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035 USA 
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I. Introduction 


The proposed National Aerospace Plane will be powered by an air-breathing, 
hydrogen-burning, supersonic ramjet (SCRAMjet) engine. The supersonic flow in 
the combustor and nozzle regions of the SCRAMjet will allow only very short res- 
idence times for the reactants, and consequently, accurate knowledge of the finite 
rates of reaction is critical to design of the engine. The rates of recombination 
processes are of special interest since these reactions can lead to considerable heat 
evolution and their rates, under conditions of high temperature and supersonic flow, 
are not well-known experimentally. One such reaction which is expected to be im- 
portant under these conditions is the recombination of two hydrogen atoms in the 
presence of N 2 : 


H + H + N 2 — * H 2 “h A^2 (1) 

The various mechanisms proposed to describe three-body recombination pro- 
cesses involve the formation of a metastable species. For reaction (1), the expected 
metastable species is either Hj (vibrationally-rotationally excited H 2 ), which can 
be stabilized by subsequent collisions with N 2 , or an HN 2 species, which can then 
react -with a hydrogen atom to yield bound H 2 : 

H + EN 2 -+ H 2 + N 2 (2) 

Calculations characterizing reaction (2) will be the subject of a future publication. 
We simply note here that these calculations indicate there is no barrier to reaction 
(2). However, in addition to the region of the potential energy surface which controls 
the abstraction process, there are three additional regions of the surface which lead 
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to stable H 2 N 2 species (cis HNNH, trans HNNH, and H 2 NN). While the HN 2 species 
has not been observed experimentally, indirect evidence for its existence has been 
obtained from experimental studies of the reaction: 

NH 2 + NO — > products (3) 

This reaction constitutes one possible step in the complicated mechanisms which 
describe the fundamental chemistry of nitric oxide both in the troposphere 1 and 
in combustion 2-5 processes. As a result, it has been studied extensively with a 
variety of experimental techniques over a range of temperatures and pressures. Early 
studies 6-15 suggested that the major pathway is the reaction: 

NH 2 + N0^N 2 + H 2 0 (4) 

Vibrationally excited H 2 0 and small amounts of the NH 2 NO complex were observed 
along with a negative temperature dependence of the rate coefficient. Gehring et 
al. 11 did entertain the possibility of an alternate pathway: 

NH 2 + NO -► N 2 + OH + H (5) 

However, a careful search using ESR techniques failed to identify any measurable 
quantity of hydrogen atoms. 

In subsequent studies 16-29 , investigators considered a reaction mechanism which 
includes three possible product channels: 

NH 2 + N0 -» N 2 + H 2 0 (6a) 
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NH 2 + NO -» N 2 + OH + H 


(6b) 


NH 2 + JVO HN 2 + OH (6c) 

Each of these experiments attempted to determine the rate coefficient of the NH 2 + 
NO reaction as well as the major pathway of this reaction. The results of the studies 
are far from consistent. Stief et al. 18 did not observe OH radicals, suggesting that 
the branching ratio for the formation of OH radicals is substantially less than 0.22 
at 298 K. Further, both ESR techniques and Lyman-alpha fluorescence could not 
detect the production of hydrogen atoms. Hall et al. 28 , using infrared absorption 
at 298 K, reported that the OH channels account for only 13 ± 2% of the reaction, 
while Dolson 29 attributed 15% of the reaction to OH radical production at 300 K. 
In contrast, Silver and Kolb 17 determined that approximately 40% of the NH 2 + 
NO reaction produced OH radicals over a temperature range of 294 - 1215 K, 
although, at room temperature, they did not observe hydrogen atoms. Andresen 
et al. 20 , operating an isothermal flow reactor at temperatures ranging from 300 K 
to 1150 K, reported that the OH branching ratio is greater than 0.65 but that 
hydrogen atoms were not detected for temperatures from 300 K to 900 K. Finally, 
Kimball-Linne and Hanson 27 observed a branching ratio for OH radical formation 
of 0.48 at 1050 K and a ratio of greater than 0.80 at 1400 K. 

Despite the wide variation in the reported branching ratios for OH radical forma- 
tion, the experimental results consistently indicate, over a range of temperatures, 
that hydrogen atoms are not produced. In light of equation (6), the observation 
of OH radicals without the production of hydrogen atoms indirectly argues for the 
existence of the rather exotic species, HN 2 . However, none of the experimental 
studies have directly detected the presence of HN 2 . Additional arguments for the 
existence of the HN 2 species have been made by researchers constructing kinetic 
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models of combustion processes 3 ’ 4 . These models include several reactions of HN 2 
as critical steps and require that this species possess a lifetime which is long enough 
to allow its participation in these reactions. 

There have been a number of previous ab initio studies 30-35 of the reaction: 

HN 2 -> H + N 2 (7) 

characterizing the geometry and energetics of the reactant, transition state, and 
products. The most recent of these earlier studies, the work of Curtiss et al. 35 , 
concisely summarizes the methodologies and basis sets used in these investigations. 
At the SCF level of theory, the dissociation energy of HN 2 to H + N 2 was deter- 
mined to have a value varying from -9.6 kcal mol -1 to +13.6 kcal mol -1 with a 
barrier to dissociation (before zero-point energy corrections) of 22.6 kcal mol -1 . 
The first two studies which treated electron correlation failed to resolve these early 
discrepancies. Vasudevan et al. 31 estimated the dissociation energy of HN 2 to H + 
N 2 at 3.8 - 4.4 kcal mol -1 with an barrier height of 10 - 15 kcal mol -1 . Casewit 
and Goddard 34 concluded that the dissociation energy is —20.0 kcal mol -1 but did 
not reinvestigate the barrier height. In their study, Curtiss et al. 35 utilized third- 
order and fourth-order Mpller-Plesset perturbation theory (MP3/MP4) to account 
for electron correlation effects and a double zeta basis set which included polariza- 
tion functions (6-31 G**) to compute a dissociation energy of —14.4 kcal mol -1 and 
a barrier of 10.5 kcal mol -1 . Using their calculated geometries and energetics in 
conjunction with Conventional Transition State Theory (including an estimate of 
quantal tunneling), Curtiss et al. then determined an approximate lifetime for the 
HN 2 radical. 

While M0ller-Plesset perturbation theory techniques work well in regions which 
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are well-described by single configuration SCF wavefunctions, difficulties do arise 
when bonds are stretched and when multiply-bonded systems such as N 2 are in- 
volved. Given the possible importance of the HN 2 species in a number of chem- 
ical systems, we decided to carry out an extensive calculation of the potential 
energy surface for the addition of hydrogen atom to N 2 , using Complete Active 
Space/externally Contracted Cl wavefunctions (CASSCF/CCI) 36 along with good 
quality Atomic Natural Orbital (ANO) basis sets 37 . We have utilized these new 
calculations to study the lifetime of the HN 2 radical, using both Conventional Tran- 
sition State Theory 38 and the formalism of W. H. Miller 39 which utilizes an Eckart 
Barrier to estimate one-dimensional quantum mechanical tunneling. 

The ab initio computational method is discussed in Section II. with the results of 
these calculations presented in Section III. Section IV. then discusses the dynamical 
methods used in this study along with an estimate of the lifetime of the HN 2 species. 
Finally, Section V. summarizes our conclusions. 

II. Ab Initio Computational Details 

Using a localized orbital description, the ground state of HN 2 is represented by 
equation (8): 


In equation (8), one NN ir bond has been broken and an NH a bond has been formed. 
Since the two bonds (one formed and one broken) are of about equal strength, a 
significant barrier is expected to occur in the entrance channel. Consequently, the 
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potential energy surface for hydrogen atom addition to N 2 has two stationary points, 
a saddle point for Rnh « 2.7 a 0 and a minimum for Rnh ~ 2.0 ao. 

ANO basis sets 37 are used in the present study. Since these basis sets are 
optimal for describing the atomic correlation effects, they have very small basis 
set superposition errors, but they are sufficiently flexible to be used in molec- 
ular calculations at both the SCF and Cl levels. The N and H basis sets are 
(13s8p6d4f)/[4s3p2dlf] and (8s6p4d)/[3s2pld], respectively, and are described in 
detail in reference 37. While these basis sets were expected to describe the molecu- 
lar potential energy curves accurately, additional calculations were completed with 
larger N, (13s8p6d4f)/[5s4p3d2f], and H, (8s6p4d)/[4s3p2d], basis sets, in which 
the additional contracted functions correspond to higher energy atomic natural or- 
bitals. These additional calculations were carried out at both the stationary point 
geometries determined with the smaller basis sets, and the H -|- N 2 geometry. 

Table I. describes the qualitative character of the CASSCF orbitals for both H 
+ N 2 and HN 2 . The la' and 2a' orbitals correspond to the N Is orbitals, which 
were not correlated in these calculations, while the 3a' and 4a' orbitals correspond 
to the N 2s orbitals. The 5a' and 9a' orbitals are the natural orbitals of the N 2 <r 
bond pair. At large Rnh, the 6a' orbital corresponds to the N 2 l7r u orbital, the 7a' 
orbital to the H Is orbital, and the 8a' orbital to the N 2 l7r g orbital. At shorter 
R-NH) there is considerable mixing of the orbitals which are appropriate in the limit 
of large Rnh- Near the HN 2 minimum, the 6a' orbital corresponds to the N 2 1tt u 
+ H Is, the 7a' orbital to N 2 l7r g + H Is, and the 8a' orbital to N 2 l7r g — H Is. 

Orbitals 4a' through 9a', and the la" and 2a" orbitals were active in the CASSCF 
calculations. This CASSCF treatment in six active a' and two active a" orbitals is 
denoted as a (62) CASSCF. The inclusion of the 4a' orbital in the active space was 
not conceptually necessary based on equation (8), but it was found that the larger 
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(62) active space led to better CASSCF convergence, even though the CASSCF 
wavefunction of this orbital contained no configurations with Cl coefficients greater 
than 0.05. In the CCI calculations, all but the N Is electrons (the la/ and 2a' 
orbitals) were correlated. The CCI list of reference configurations was obtained by 
selecting the unique configurations from the H + N 2 region, the saddle point region, 
and the HN 2 minimum region whose Cl coefficients in the CASSCF wavefunction 
were greater than 0.05. The multireference analogue of Davidson’s correction 40 was 
computed to estimate the importance of higher-order excitations. This correction 
to the CCI energy is given by 


Q = AE(1 - C’)/Cj (9) 

where AE is the Cl energy minus the reference energy, and Cq is the sum of the 
squares of the coefficients of the reference configurations in the Cl wavefunction. 

The calculations using the smaller basis set were carried out on the NASA Ames 
Cray X-MP/48. These calculations used the MOLECULE - SWEDEN 41 ’ 42 system 
of programs. The calculations using the larger basis set were carried out at Sandia 
National Laboratory on a Cray X-MP/416. 

III. Ab Initio Results 

Table II. lists the parameters of the N 2 molecular potential energy surface which 
were computed with the [4s3p2dlf] and the [5s4p3d2f] basis sets. In order to be 
consistent with the HN 2 analysis, the diatomic calculations included the N 2p de- 
rived orbitals and the N 2 2<7 U orbital in the active space. The CCI reference fist was 
larger than in the HN 2 case, consisting of all the configurations generated from the 
N 2p derived levels. As the Table indicates, the smaller basis set yields values for 
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R e , u>e, and Do which differ from the corresponding experimental 43 quantities by 
0.022 ao, 73 cm -1 , and 13 kcal mol -3 (based on a Dunham analysis of a quadratic 
fit near R e ). The value of D 0 computed in this work differs from that previously 
calculated with a segmented basis set of the same size 44 by 0.5 kcal mol -1 . This 
earlier study 44 concluded that the error resulting from the externally contracted Cl 
procedure is expected to be only 0.16 kcal mol -1 (determined by comparing CCI 
+ Q and SDCI + Q results). Consequently, limitations in the basis set itself are 
responsible for the remaining difference between the value of Do calculated in this 
work and the experimental value. 

The value of u> e for N 2 computed in this study is 248 cm -3 smaller than that 
obtained from the MP3/6-31G** calculations of Curtiss et al. 35 , and underesti- 
mates the experimental harmonic frequency by approximately 3.1% compared with 
an overestimate of approximately 7.4% by the M0ller-Plesset perturbation theory 
calculations. Curtiss and co-workers suggested that the discrepancy between their 
calculated N 2 frequency and the experimental value is a consequence of limitations 
in the 6-31G** basis set. This conclusion is consistent with the earlier work of 
Siegbahn 36 . Additionally, as Bartlett and Purvis 45 , and more recently, Laidig et 
al. 46 have discussed, the multiply-bonded N 2 system cannot be represented sat- 
isfactorily using only a single-reference, many-body perturbation theory (MBPT) 
technique. Consequently, it appears that limitations in both the basis set and 
the computational methodology have contributed to the anomalous N 2 frequency 
reported by Curtiss and co-workers. It should be noted, however, that even the 
CASSCF/CCI calculations reported in this work, utilizing a good quality ANO 
basis set, fail to represent precisely the experimental values of both D 0 and the 
harmonic stretching frequency of the N 2 triple bond system. 

The larger basis set, [5s4p3d2f], reduces these errors to 
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The first three columns of Table III. summarize the computed geometries and 
harmonic vibrational frequencies at the stationary points of the H -+■ N 2 potential 
energy surface (using the [4s3p2dlf/3s2pld] basis set). The values of Rnn and o> 2 
at the H + N 2 minimum were derived from supermolecule calculations and are very 
similar to the values of R e and u> e obtained for N 2 with the [4s3p2dlf] basis set (see 
Table II). These results differ because the diatomic calculations used all possible 
reference configurations from the active space to estimate electron correlation, while 
the supermolecule calculations utilized a set of selected reference configurations from 
the same active space. Energetically, the HN 2 species lies 3.01 kcal mol -1 above the 
H + N 2 asymptote, while the saddle point for hydrogen atom addition is 15.16 kcal 
mol -1 above H + N 2 . Thus, there is a barrier of 12.15 kcal mol -1 to the dissociation 
of HN 2 to H + N 2 . The addition of zero- point energy corrections reduces this barrier 
to 8.49 kcal mol -1 . Single-point calculations with the [5s4p3d2f/4s3p2d] basis set 
(at the stationary point geometries determined with the [4s3p2dlf/3s2pld] basis set) 
indicate that the HN 2 species and the saddle point are 3.73 kcal mol -1 and 15.01 
kcal mol -1 , respectively, above the H + N 2 asymptote, thus yielding a dissociation 
barrier for HN 2 of 11.28 kcal mol -1 . 

The harmonic vibrational frequencies of Table III. are based on normal mode 
analyses performed at the stationary points of the H + N 2 potential energy sur- 
face. u> 1 , o> 2 , and u> 3 are the frequencies of the normal modes which correspond 
principally to the NH stretching motion, the NN stretching motion, and the HNN 
bending motion. The value of u> 2 decreases from 2303 cm -1 at the H + N 2 asymp- 
tote, to 2072cm -1 at the H — N 2 saddle point, and finally, to 1583 cm -1 at the 
HN 2 minimum. Simultaneously, the value of (the frequency which essentially 
corresponds to the reaction coordinate for hydrogen atom addition to N 2 ), changes 
from an imaginary value (1662t cm -1 ) at the H — N 2 saddle point to a value of 2744 
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cm -1 at the HN 2 minimum. Finally, the value of increases from 771 cm -1 at 
the H — N 2 saddle point to 1070 cm -1 at the HN 2 minimum. These changes are 
consistent with the breaking of on NN 7 r bond and simultaneous formation of an 
NH a bond as the HN 2 species forms. 

Table III. also compares the results obtained in the present study at the stationary 
points of the H + N 2 surface with those reported by Curtiss et al. 35 . The current 
calculations suggest that the HN 2 species is energetically more stable than predicted 
by the earlier work, lying only 3.01 kcal mol -1 rather than previously calculated 14.4 
kcal mol -1 above the H + N 2 asymptote. In addition, the barrier to dissociation 
(before zero-point energy corrections are added) is currently found to be 12.15 kcal 
mol -1 compared with the earlier value of 10.6 kcal mol -1 . Much of this difference 
can be attributed to the larger basis set employed in the CASSCF/CCI study, since 
the basis set used by Curtiss and co-workers did not contain f functions on nitrogen 
or d functions on hydrogen. The d and f functions, which describe higher angular 
momentum states, are essential to the calculation of an accurate D e value for the 
HN 2 species. 

As Table III. indicates, there are significant differences among the normal mode 
frequencies calculated in the current work and in that of Curtiss and co-workers. 
The values of u> 3 (primarily the HNN bending motion) computed in the two studies 
are in good agreement at both the HN 2 minimum and the saddle point of the H + 
N 2 potential energy surface. However, this is not true of the other two normal mode 
frequencies. The calculated values of a>j (primarily the NH stretching motion) and 
o> 2 (primarily the NN stretching motion) differ substantially at the stationary points. 
Further, while the study of Curtiss and co-workers reports a puzzling increase in 
the u> 2 frequency at both the saddle point and at the HN 2 minimum compared to 
its value in the N 2 molecule despite a lengthening of the NN bond, the u >2 frequency 
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does not exhibit such anomalous behavior in the CASSCF/CCI calculations. In the 
current study, as the NN bond lengthens, the w 2 frequency decreases monotonically. 
The unexplained behavior of the u> 2 frequency in the work of Curtiss et al. suggests 
a dramatic breakdown of the Mpller-Plesset perturbation expansion. As noted 
above, this effect is not surprising since the N 2 triple bond requires a multireference 
description to correctly represent the elongation of the NN bond. 

Finally, we determined the minimum energy path (MEP) for the addition of H 
atom to N 2 by computing a series of two-dimensional grids defined in terms of Rnn 
and 9 (the HNN angle) for various values of Rnh (Rnh = 4.5, 4.0, 3.5, 3.0, 2.8, 2.6, 
2.4, 2.2, 2.0, and 1.8 ao). A quadratic polynomial in Rnn and 6 was then fitted 
to each of these calculated grids (using a linear least-squares procedure), and the 
minima of the resulting polynomials located. While all the computed energy values 
are contained in the appendix, Table IV. lists the coordinates, the corresponding 
potential energy values, and the normal mode frequencies orthogonal to the MEP 
(computed from the quadratic polynomials). Figure 1 depicts the energy along the 
MEP as a function of the Rnh distance. At large Rnh distances, the H atom initially 
approaches N 2 at an angle of rs 119°, with 6 decreasing to ^ 115° as Rnh assumes 
its equilibrium value in HN 2 . Simultaneously, the NN bond gradually elongates 
as the Rnh distance decreases. These features of the MEP are similar to those 
observed in the addition of H atom to 0 2 47 . 

TV. Dynamical Results 

From the point of view of reaction dynamics, the lifetime of a species is perhaps 
the most elementary molecular property open to investigation. In the case of an 
extremely reactive intermediate, this parameter plays a crucial role in determining 
the ability of the species to participate in subsequent chemical reactions. We now 
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present a relatively straightforward, two-step analysis of the HN 2 lifetime, utilizing 
both Conventional Transition State Theory in its thermodynamic formulation, and 
an Eckart barrier formalism which estimates the role of one-dimensional quantum 
mechanical tunneling in the unimolecular decay process. This analysis uses infor- 
mation from the stationary points and the MEP of the molecular potential energy 
surface, but does not attempt to describe the microscopic details of the HN 2 decom- 
position process. An investigation of these details requires a more comprehensive 
representation of the potential energy surface and is left to future work. Rather, 
our intent is to obtain a new estimate of the HN 2 lifetime based on the current 
CASSCF/CCI calculations, to compare this estimate with a previous determina- 
tion of the HN 2 lifetime, and to assess the implications of these lifetime estimates 
for the chemistry of the HN 2 species. 

The upper portion of Table V. contains a summary of the thermodynamic proper- 
ties of the H + N 2 — * HN 2 reaction based on information determined at the station- 
ary points of the potential energy surface and calculated with statistical mechanics. 
This data is divided into two sections, the first reflects the changes in the thermo- 
dynamic state functions between the HN 2 minimum and the H + N 2 asymptote, 
while the second tabulates the same changes between the HN 2 minimum and the 
transition state (quantities marked by *). In both cases the HN 2 minimum has been 
chosen as the zero of energy. The data is from three distinct sources: (1) column 
one is based on the current CASSCF/CCI calculations using the [4s3p2dlf] basis 
set; (2) column two contains results from recent Bond Additivity Correction with 
fourth-order Mpller-Plesset perturbation theory (BAC-MP4) calculations of Melius 
and Binkley 48 ; (3) the remaining columns summarize the results of the MP3/MP4 
calculations of Curtiss and co- workers 35 . 

Before discussing the HN 2 lifetimes computed from the thermodynamic data us- 
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ing Conventional Transition State Theory, several characteristics of the data itself 
deserve comment. While columns three and four list the results reported by Curtiss 
et al. (Table 6. of reference 35.), columns five and six contain the same quantities 
calculated directly from the energetics, geometries, and frequencies of reference 35. 
In the first section of the data (relating changes between the HN 2 minimum and 
the H + N 2 asymptote), there are minor discrepancies (several tenths of a kcal 
mol -1 ) between the reported and the calculated values of AH at both 0 K and 
298.15 K. These differences arise from the apparent use by Curtiss and co-workers 
of both an N 2 frequency and bond length near the experimental values (2359 cm -1 
and 2.075 ao) to prepare their reported data, compared with the MP3/6-31G** 
harmonic frequency of 2534 cm -1 and computed bond length of 2.090 ao used 
to generate the corresponding calculated MP3/MP4 values. We have utilized the 
MP3/6-31G** harmonic frequency and bond length in preparing the calculated re- 
sults to allow a more direct comparison with the thermodynamic quantities based 
on the CASSCF/CCI calculations. There is, however, a more serious discrepancy 
between the reported and calculated values of AS; they differ in both magnitude 
and sign. While two different computational approaches have yielded consistent 
MP3/MP4 values of AS which agree with both the CASSCF/CCI and BAC-MP4 
(computed independently) results, we have been unable to identify a reason for the 
differences observed in Table V. Since the a value of AS is not used in the Transi- 
tion State Theory calculations of the HN 2 lifetime, the conclusions of Curtiss and 
co-workers are not affected by this apparent anomaly in the reported AS value. 
The remaining thermodynamic data of Table V. (relating changes between the HN 2 
minimum and the transition state) exhibit some variations between the reported 
and calculated MP3/MP4 values which do affect the computed HN 2 lifetimes, but 
only in a relatively minor way. Again, the conclusions of Curtiss et al. are not 
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affected by these small changes. 

The lower portion of Table V. summarizes the classical (that is, no quantum 
mechanical tunneling) HN 2 lifetimes computed with conventional Transition State 
Theory from the expression: 



where h is Planck’s constant, k is Boltzmann’s constant, R is the gas constant, and 
T is temperature (equation (10) is simply the reciprocal of the Transition State 
Theory rate coefficient for the unimolecular decay of HN 2 ). The values reported by 
Curtiss and co-workers are consistent with the calculated MP4 results listed in the 
Table, while the lifetimes predicted by the CASSCF /CCI data are approximately 
two orders of magnitude greater at 298.15 K and approximately 19 orders of magni- 
tude greater at 30.0 K. These differences are a simple consequence of the disparate 
barriers to unimolecular decay predicted by the two sets of ab initio data. At this 
approximate level of theory it is inappropriate to attach significance to these dif- 
ferences or to the precise numerical values, and the general conclusions of the two 
studies, based on conventional Transition State Theory calculations, are the same. 
First, the lifetime of the HN 2 species at room temperature is very short (on the or- 
der of tenths to thousandths of a microsecond); secondly, at very low temperature, 
the predicted very long lifetime of H1V 2 may permit experimental observation of the 
species. 

Up to this point in our analysis, the HN 2 species has been treated in a purely 
classical manner, ignoring quantum mechanical effects. This approach has allowed a 
direct comparison with a portion of the study by Curtiss and co-workers. However, 
since the ab initio calculations discussed earlier indicate that the HN 2 minimum 
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lies at least 3 kcal mol -1 above the H -f N2 asymptote, every rovibrational state 
supported by the HN2 well is only quasi-bound. Consequently, quantum mechanical 
tunneling plays an essential role in the unimolecular decay of HN2. Curtiss and co- 
workers accounted for these quantal effects by computing the tunneling frequency 
of a bound particle through a parabolic barrier. In this study we have used an 
Eckart function to model the vibrationally adiabatic barrier to HN2 dissociation, 
computing HN2 lifetimes based on an estimate of the one-dimensional quantal tun- 
neling through that barrier. The vibrationally adiabatic potential was determined 
by adding the zero-point energy contribution of the NN stretching motion and the 
HNN bending motion (frequencies listed in Table IV.) to the computed potential 
energy along the MEP. Further, because all of the rovibrational states of the HN 2 
species appear to be only quasi-bound, it does not seem appropriate to compute uni- 
molecular lifetimes by assuming a thermally- averaged Boltzmann distribution (as 
does the thermodynamic formulation of Transition State Theory) created by mul- 
tiple collisions of the HN2 species with other molecules in the gas-phase. Instead, 
the formalism of W. H. Miller 39 was utilized to calculate unimolecular lifetimes as 
a function of the energy of each HN2 quasi-bound vibrational state (microcanonical 
lifetimes). For the purposes of this study, we have assumed that J = 0 (that is, the 
HN 2 species is rotationally cold), although the formalism is not limited to this case 
alone. 

Table VI. summarizes the calculated unimolecular lifetimes of HN 2 as a func- 
tion of the energy of each quasi-bound vibrational state. In order to provide a 
careful estimate of the quantal tunneling effects, the calculations reported in Table 
VI. represent the results of three distinct Eckart Barrier models (designated Eckart 
Barrier- 1 , Eckart Barrier- 2 , and Eckart Barrier- 3 ). In these models, the three pa- 
rameters which completely determine the Eckart function, the height of the barrier 
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with respect to HN 2 (V 0 ), the height of the barrier with respect to H + N 2 (Vi ), and 
the imaginary frequency (which determines the curvature at the top of the barrier) 
(u>*), were systematically varied in an attempt to reflect accurately the properties 
of the HN 2 potential energy surface. 

In the first model, Eckart Barrier-1, no effort was made to fit numerically the 
Eckart function to the shape of the vibrationally adiabatic potential. Rather, the 
parameters Vo, Vj, and u>* were set equal to the zero- point-corrected (adding the 
zero-point energy of the NN stretch and the HNN bend) ab initio potential en- 
ergy values and to the calculated imaginary frequency. Table VI. reports three 
calculations of the HN 2 lifetime as a function of energy using the Eckart Barrier- 1 
model: (1) a calculation using the ab initio results of the [4s3p2dlf] basis set, with 
V 0 = 12.417 kcal mol -1 , Vi = 15.931 kcal mol -1 , and u > * = 1662.2t cm -1 ; (2) a 
calculation using the ab initio results of the [5s4p3d2f] basis set (at the [4s3p2dlf] 
geometry), with Vo = 11.549 kcal mol -1 , Vi = 15.777 kcal mol -1 , and u* = 1662.2t 
cm -1 ; (3) a calculation using the MP4/6-31G** results of Curtiss and co-workers 
with Vo = 10.310 kcal mol -1 , Vj = 26.350 kcal mol -1 , and w* = 1533i cm -1 . Fig- 
ure 2. compares the Eckart Barrier- 1 model to the vibrationally adiabatic potential 
computed with the [4s3p2dlf] basis set. As Ihble VI. indicates, the MP4/6-31G"* 
potential energy surface of Curtiss and co-workers supports only two quasi-bound 
vibrational states, while the potential energy surfaces computed in this work with 
the [4s3p2dlf] basis set and the [5s4p3d2f] basis set support six and five quasi- 
bound vibrational states, respectively. The Table also demonstrates the dynamical 
effect of using a larger ab initio basis set on the HN 2 lifetimes which range from 
6.37 x 10 -e seconds to 1.76 x 10 -13 seconds for the six vibrational levels of the 
[4s3p2dlf] basis, and from 1.25 x 10 -9 seconds to 1.02 x 10 -13 seconds for the five 
vibrational levels of the [5s4p3d2f] basis. Finally, the lifetime computed for HN 2 
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dissociation from the lowest vibrational state of the [4s3p2dlf] potential energy sur- 
face is approximately two orders of magnitude greater than the lifetime calculated 
for the unimolecular decay of HN 2 from the lowest vibrational state supported by 
the MP4/6-31G** potential energy surface. 

However, as Figure 2. indicates, the Eckart Barrier-1 model both significantly 
overestimates the width of the barrier and does not closely approximate the shape 
of the vibrationally adiabatic potential. We therefore attempted to define a single 
model which is able to represent accurately both the width of the barrier and 
the shape of the vibrationally adiabatic potential. This proved to be impossible. 
Consequently, two additional models were constructed by numerically fitting (in a 
non-linear, least-squares sense) an Eckart function to the vibrationally adiabatic 
potential determined with the [4s3p2dlf] basis set. In the first of these new models, 
designated Eckart Barrier-2, Vo was held fixed at the value used in the Eckart 
Barrier-1 model, while Vj and u >* were treated as adjustable parameters. The 
fitted Eckart Barrier- 2 model is characterized by the parameters V 0 = 12.417 kcal 
mol -1 , Vj = 12.997 kcal mol -1 , and u>* = 2441. 3i cm -1 , but, as Figure 3. indicates, 
it significantly underestimates the width of the barrier. This model does, however, 
possess the attribute of closely approximating the shape of the HN 2 side of the 
vibrationally adiabatic potential. The second of the new models, identified as Eckart 
Barrier-3, was constructed by holding both Vo and Vj fixed at the values which 
characterize the Eckart Barrier-2 model (hence, Vo is the same value in all three 
models considered in this work) and treating as the only adjustable parameter 
in the non-linear least-squares fitting procedure. The resulting Eckart Barrier-3 
model is defined by the parameters Vo = 12.417 kcal mol -1 , Vj = 12.997 kcal 
mol -1 , and u>* = 1715. 6i cm -1 , and is compared in Figure 4. to the vibrationally 
adiabatic potential. While the Eckart Barrier-3 model overestimates the width of 
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the barrier, it does succeed in closely approximating the shape of H + N 2 side of 
the vibrationally adiabatic potential. 

The two models provide a lower bound (Eckart Barrier-2) and an upper bound 
(Eckart Barrier-3) for the HN 2 unimolecular lifetime of each quasi-bound vibrational 
state associated with the [4s3p2dlf] basis set. These bounds are summarized in the 
fined two columns of Table VI. Note that as the top of the barrier is approached, 
the widths of the two model barriers become increasingly similar. Consequently, 
the difference between the two lifetimes computed for each quasi-bound vibrational 
state decreases with increasing energy. In the case of the lowest-lying quasi-bound 
vibrational state, the two computed lifetimes differ by approximately a factor of 66; 
for the highest-lying state, they differ only by approximately a factor of 1.4. The 
HN 2 lifetimes predicted by these two models, which include the effects of quantum 
mechanical tunneling, suggest that the HN 2 species is extremely short-lived, even 
at very low temperatures. As a result, we expect that it will be difficult to observe 
the HN 2 species experimentally, and that the role of HN 2 in both termolecular 
recombination processes and combustion processes will be dramatically limited by 
its very short lifetime. 

V. Conclusions 

This study has focused on two aspects of the addition of atomic hydrogen to N 2 : 
first, a characterization of the minimum energy path for the addition reaction, and 
secondly, a determination of the unimolecular lifetime of the HN 2 species. While 
this work has presented neither a global description of the potential energy surface 
for the addition reaction nor a discussion of the microscopic details of the dynamics, 
it has reached several important conclusions: 

(1) The HN 2 minimum is found to lie approximately 3.01 kcal mol -1 above the 
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H 4- N 2 asymptote, compared to 14.4 kcal mol -1 calculated earlier by Curtiss 
and co-workers. 

(2) While this earlier study computed a zero-point-corrected barrier to dissociation 
of approximately 5.8 kcal mol -1 , the current work has estimated this barrier 
to be approximately 8.5 kcal mol -1 . 

(3) The current calculations indicate that six quasi-bound vibrational levels exist 
in the HN 2 well, compared to the two quasi-bound levels predicted in the earlier 
study of Curtiss and co-workers. 

(4) The current work has used both a larger basis set and an improved methodology 
to calculate a new estimate of the NN harmonic stretching frequency in the N 2 
and HN 2 species. 

(5) We estimate that the lifetime of the lowest-lying quasi-bound vibrational state 
of HN 2 to lie between 8.8 x 10 -11 seconds and 5.8 x 10 -9 seconds. This 
predicted lifetime is between a factor of two and a factor of one hundred longer 
them the lifetime of the HN 2 species predicted by Curtiss and co-workers. 

(6) However, we conclude, in concert with the conclusions of Curtiss and co- 
workers, that it will be very difficult to observe the HN 2 species experimen- 
tally, and that the predicted short lifetime of the HN 2 species will severely 
circumscribe its role in both termolecular reaction mechanisms and proposed 
combustion mechanisms. 
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Table I. Qualitative Character of the CASSCF Orbitals for H + N 2 — ► HN^ 
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Table II. Computed Potential Curve Parameters for N 2 . 
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Table III. Stationary Points on the II T N 2 Surface.' 
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Table IV. Computed Minimum Energy Path for II + N 2 .' 
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Table V. Thermodynamic Properties and Unimolecular Lifetimes for IIN 
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Table VI. Calculated Unimolecular Lifetimes" of IIN 
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FIGURE CAPTIONS 


Figure 1.: Interaction energy for the addition of a hydrogen atom to N 2 as a 
function of th along a minimum energy path (MEP) obtained by optimizing Rnn 
and 9 at each Rnh value. The energies are from CASSCF/CCI -f Q calculations 
(see text). 

Figure 2.: Eckart Barrier-1 model. The bold solid curve is the vibrationally adi- 
abatic barrier calculated from the CASSCF/CCI + Q data along the minimum 
energy path (MEP). The solid curve is an Eckart function. The horizontal dashed 
curves represent the quasi-bound vibrational states of the HN 2 well. In the separable 
harmonic approximation, the energy of each vibrational state is given by E( ni)niin3 ) 
= (ni + |)ht/j. Each state is identified by a set of quantum numbers (n 3 , n 2 , 

n 3 ), where n 3 = HN stretch quantum number, n 2 = NN stretch quantum number, 
n 3 = HNN bend quantum number. 

Figure 3.: Eckaxt Barrier-2 model. The bold solid curve is the vibrationally adi- 
abatic barrier calculated from the CASSCF/CCI + Q data along the minimum 
energy path (MEP). The solid curve is an Eckart function. The horizontal dashed 
curves represent the quasi-bound vibrational states of the HN 2 well. In the separable 
harmonic approximation, the energy of each vibrational state is given by E( niin , in3 ) 
= (ni + |)hi/i. Each state is identified by a set of quantum numbers (n 3 , n 2 , 

n 3 ), where n 3 = HN stretch quantum number, n 2 = NN stretch quantum number, 
n 3 = HNN bend quantum number. 

Figure 4.: Eckart Barrier-3 model. The bold solid curve is the vibrationally adi- 
abatic barrier calculated from the CASSCF/CCI + Q data along the minimum 
energy path (MEP). The solid curve is an Eckart function. The horizontal dashed 
curves represent the quasi-bound vibrational states of the HN 2 well. In the separable 
harmonic approximation, the energy of each vibrational state is given by E(„ liIl2 in3 ) 
= (ni + |)hi/i. Each state is identified by a set of quantum numbers (n 3 , n 2 , 

n 3 ), where n 3 = HN stretch quantum number, n 2 = NN stretch quantum number, 
n 3 = HNN bend quantum number. 




3.5 4.0 4.5 5.0 

BOHR 





BOHR 




in 

in 



1 


10 n 1VDM 


BOHR 




BOHR 





